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Abstract 

6\". 

We describe the finite volume effects of CP-odd quantities, such as the neutron 
electric dipole moment and the anapole moment in the ^-vacuum, under different topo- 
logical sectors. We evaluate the three-point Green's functions for the electromagnetic 
current in a fixed non-trivial topological sector in order to extract these CP-odd ob- 
servables. We discuss the role of zero modes in the CP-odd Green's function and show 
that, in the quenched approximation, there is a power divergence in the quark mass 
for CP-odd quantities at finite volume. 
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CP-violation is still not well understood. The baryon asymmetry in the universe suggests 



that there is another source of CP-violation besides the CP phase in the CKM matrix as 
measured from the K and B meson decays.. The recent discovery of non-zero neutrino mass 
has prompted suggestions of CP-violation in the leptonic sector. On the other hand, it 
is known for many years that QCD has a 9 term which gives rise to strong CP-violation. 
However, its consequence in quantities such as the neutron electric dipole moment has not 
been discovered. This puts the limit of 9 to be less than 10 -9 . Next generation of experiments 
will further push down the limit and hopefully will detect it at some stage. The relation 
between the neutron electric dipole moment due to the 9 term can be calculated via lattice 
QCD and there are several attempts to do so [U EJ [3j Hj. Since 9 in QCD is small, one 
can compute the CP-odd part of the Green's function as the first derivative w.r.t. 9, i.e. 
i9(QO)g =0 = i9^2 1/ P v u(0) v . This requires a weighted sum of the Green's function (0) v in 
different topological sectors. In view of the fact that practical Monte Carlo simulations can 
get stuck in a particular topological sector, how is one to evaluate CP-odd quantities such as 
the neutron electric dipole moment and anapole moment in this case? The situation of being 
stuck in one topological sector is getting more serious as the gauge field becomes smooth, 
the quark mass is approaching the physically small u/d mass and the continuum limit is 
approached. Since simulating overlap fermion with HMC faces the difficulty of tackling 
discontinuity of the sign function associated with topological change, it has prompted JLQCD 
to add a extra fermion action to prevent the topological change to avoid this numerical 
difficulty [5j. Recently it is shown that gauge action can be obtained from the overlap 
operator as TrD ov J6j [7] and can be simulated together with the overlap fermion [8j. The 
first attempt in implementing it in HMC with the rational approximation has revealed that 
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even for heavy quark (ma = 0.5) for a range of lattice spacing a = 0.06 — 0.15 fm, no 
topological change has been encountered [9]. These initial MC runs have raised the issue of 
ergodicity. Nevertheless, if the Monte Carlo algorithm is ergodic, all the physical quantities 
including the CP-odd quantities can be obtained within a fixed topological sector provided 
finite volume correction is taken into account. This has been formulated and studied for 
CP-even quantities such as the mass [10] and it has been recently extended to CP-odd 
quantities and shown that topological susceptibility can be obtained in a fixed topological 
sector [11]. In the present work, we shall concentrate on the three-point function and discuss 
the electromagnetic form factors and the neutron electric dipole form factors in the fixed 
topological sector. In the above expression to extract NEDM, the Green's function (0) v is 
weighted with v. Since the average of topological charge \v\ grows with yV, in what way 
does the above expression converge at large volume? We shall address this point. We shall 
point out that zero modes are part of the physical spectrum and they contribute to the 
leading 1/V behavior in the CP-odd case. As such, the quenched calculation of NEDM and, 
for that matter, other CP-odd matrix elements suffer power divergence in quark mass at the 
chiral limit. Finally, we will discuss the possible relation between the anomalous magnetic 
form factor and NEDM of the nucleon. 

2 6 Vacuum and Topological Sectors 

We shall start with a brief review of the relation between the 9 vacuum and the associated 
topological sectors in QCD. The QCD partition function in a 9 vacuum on a torus in the 
continuum formulation is written as 

Z{9) = J VAV^V^e- s[A ^ ] e ieQ[A \ (1) 
where the topological charge operator Q[A] is defined as 

Q[A] = vhJ d * xTr ( G ^»u). ( 2 ) 

For a torus of V = L 3 L t with periodic or anti-periodic boundary condition, the topolog- 
ical charge operator defined with Gisparg- Wilson fermion, such as the overlap fermion [13] 
will have integer charges v eigenvalues on smooth gauge configurations. Since the gauge 
configurations occurring in the path-integral in Eq. ([1]) are characterized by an integer, the 
partition function of the 9 vacuum is a sum of partition functions characterized by this 
integral topological charge v 

u=oo 

Z{9)= z l6v Z v (3) 

v=— oo 

With v being an integer, Z(9) is periodic, i.e. Z{9 + 2tc) = Z(9), and therefore Z v is the 
Fourier transform of Ziff) in the range 9 G [— 7r, tt] 

Z v = ±- [ n deZ{6)e- iQv (4) 
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The partition function Ziff) in Eq. ([I]) can be expressed in terms of the energy eigenstates 



Z(6) =Y.e~ Ven{8 \ (5) 

n 

where e n (6) is the energy density of the nth eigenstate. When the spatial volume is large, the 
energy of the 9- vacuum is Eq = Veo(9). For 9 = tt, the derivative of Z{9) is discontinuous at 
9 = ±7i and CP symmetry is spontaneously broken [H]. For small 9, we assume that Z{9) 
is analytic and, on account of the fact that Z(9) = Z(—9), Cq{9) can be expanded around 
9 = " 

e {9) = \ X te 2 + i9 4 + ... (6) 

where 

Xt = ~y^, (7) 

is the topological susceptibility. At low temperature with L t (Ei(9) — E (9)) ^> 1 (L t = l/T), 
the vacuum state dominate and, therefor, 

Z{9) = e - Veo{e) = e~^ e \l - 7 fl 4 + 0(9 6 )). (8) 

When the temperature is not low compared to the energy gap which is the pion mass, 
the partition function is the Boltzmann factor which contains contributions from the excited 
states depending on the temperature. In this case, Eqs. ([8]) and © still hold provided 
that the three- volume is large and the vacuum energy density is replaced by the free energy 
density E(T, 9) (151 [P2] . We shall only consider the low temperature situation in this work. 

At low temperature and small topological charge with \ v\/yV fixed, the partition function 
for the topological sector v in Eq. (jlj) becomes 

Z v = — r d9e- VeoW e~ iev . (9) 

where the fixed v partition function Z v is dominated by the 9 vacuum with smaller energy 
density 60(6*). From Eq. ([H]), it is clear that 60(6*) has the lowest energy density at 9 = 0. 
Thus, at large V and keeping \v\/yV fixed, one employs a saddle point expansion around 
9 = to obtain 

„ 1 v' 2 , 37 67Z/ 2 7Z/ 4 , . . 

27iVxt Vxt V 2 xt V s xt 

for a small 7. 

It has been shown that physical quantities in the 9 vacuum, such as hadron masses, can 
be extracted from fixed topological sectors [TD]. Consider the n-point Green's function with 
operators Oi in the 9 vacuum 

G{6) = (0 1 2 ...O n )(9) 

= / VAV^Vm0 1 2 ...O n e~ s[A ^'^e teQ[A \ (11) 

Z [9) J 
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The corresponding Green's function in the topological sector v is 

G v = {O x 2 ...O n ) v = ^r^- T d9Z(9)G(9) e~ i0v . (12) 

Using the same saddle point approximation that was applied to Z v in Eq. ffTUl) to G u , 
Eq. f|T2|) becomes 



Gu = G{d s ) - p^e s G\e s ) + -L-G»(e.)(i - + ^- 2 G'"\e s ) + e>(-U (13) 

Vxt ZVxt Vxt 8V 2 x 2 V s 

where the pure imaginary 9 at the saddle point is 

This has been applied to two-point correlation function to show that the mass which has 
the form M(9) = M(0) + |M"(0)# 2 in the 9 vacuum is calculated in topological sector v 
is [Ha to be 

lM"(0) h v 2 , 12 39z/ 2 14z/ 4 z/ 6 

Here, we have taken the full volume dependence into account. From this expression, we see 
that the mass of the ^-vacuum M(0) can be obtained from M v at several volumes or with 
several topological sectors, or the combination thereof. Similarly, the topological suscepti- 
bility can be extracted from Eq. (j!5p with several v sectors. It is easy to prove that, upon 
summing over the topological sectors with the probability Z v in Eq. (TTOT) . the mass in the 
9 = vacuum is recovered, i.e. 

£>„M„ = M(0). (16) 



J2u Z v 



M v in Eq. ( |T5l) shows that the physical quantity measured in the fixed topological sector 
is the same as that in the 9 = vacuum with finite volume corrections. When the volume 
approaches infinity, the difference goes away, so does the difference when results from dif- 
ferent topological sectors are added up and weighted by the probability of occurrence of the 
topological sector. 

In the present work, we extend this analysis to three-point functions with the electro- 
magnetic form factor in order to extract the neutron electric dipole moment and the anapole 
moment. For the case that \9 S \ is small, we can expand G v around 9 = and obtain 

G ~ G(0) + (- lU )G'(0)(l 12l ) + G " {0) \l ^ 127 1 

-ivG>»(0) g^O) 1 
+ 2{V X tY + 8V*tf +U{ V* } - [U) 
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When the Green's function is CP-even(odd), G(9) is even(odd) in 8, the above expression 
holds separately for the even-odd cases as 

Qeven _ G(Q) , G»(0) * 12 7 G'"'(0) 1 



= ^(o)(i-^)+^|+o(i). 



It is now easy to see why the expression for evaluating the CP-odd quantities that was 
alluded to in the introduction has the correct large V behavior. To evaluate CP-odd quan- 
tity at small 6, it is sufficient to consider the derivative of the corresponding path-integral 
representation of the Green's function in Eq. (fTTj) w.r.t. 6, i.e. 



G'(0) = i(QO) = — — £ Z u v(0)u. (20) 

Plugging Z v from Eq. ffTUj) and {0) v from Eq. ffT9l) into Eq. fl20l) . one recovers G'(0) to order 
^. It is straight-forward to show that had all the V terms in Eq. f|T9|) been kept, one would 
have recovered G'(0) exactly as is for M(0) in Eq. (fT6l) . Since the leading 1/V contribution 
of (0) v is 2^G'(0) in Eq. pi]), its contribution to G"(0) in Eq. (j20j) is oc Tn,Z v ^ which 

is independent of V as it should. In other words, the apparent y/V dependence from the 
explicit v in Eq. ([20]) does not lead to \fV divergence. Similarly, it is easy to check that the 
Vxt and v dependence associated with G"(0) and G"'(0) in Eqs. ( |T8i) and ( |T9i) are correct 
to reproduce the the second and third derivatives of the path-integral formulation of G{9). 



3 Neutron Electric Dipole Moment 

We want to extend the study of two-point function to three-point function with the electro- 
magnetic form factors in the nucleon in order to extract the neutron eclectic dipole moment 
and anapole moment. The form factor for the electromagnetic current J™ = i J2 q ^q^ql^q 
in the nucleon is defined as 

ti>', s'\J* M \p, s) = u s ,{p')TM 2 >s{p), (21) 

where 

q F^(q 2 ) 
T^(g 2 ) = il^F^q 2 ) - ia llu —^F 2 (q 2 ) + {il^q 2 - 2m N ^q^)F A (q 2 ) - a^q^ , (22) 

with F^{q 2 ) /2rriN and Fa being the electric dipole and anapole form factors and = 
^[jn^Jv]- The electric dipole moment is 



4, = (23) 
2m N 
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Note that due to the charge conservation of the electromagnetic current d^J^ = 0, one has 
the relation 

(p', s'\d^ M \ P , s) = uAp')q^,(q>s( P ) = o. (24) 

The electromagnetic form factors, the neutron electric dipole and anapole form factors of 
the nucleon can be obtained from the three-point Green's function [T6l [2] . We shall extend 
it to the case at fixed topology. Consider the three-point Green's function in a 9 vacuum 

G» JN (q,t,t f ) = (XN(p',t f )J* M (q,t)x N (p,0)) e , (25) 

where q = p' — p and xn is the nucleon interpolation field. At large time separation, i.e. 
tf — t ^> 1 and f > 1, the Green's function is dominated by the lowest state which is the 
nucleon 



E ' N e (*/-*) e ~E N et 

x ^^Ixivl^Cp', S 0><^(^, sOl^^l^Cp, 5)><A^(p, 5)1x^1^) (26) 



G N jn(0, q,t,tf) ^^e-V^^e^ 



where \9) = e t6Q \0) and E' NB = \jp' + m 2 NB , E N e = ^Jp* + m 2 N g. The matrix elements for 
the interpolation fields are 

(9\ X n\N (p,s)) = Z e N u%(p,s) 
(Ne(p', s)\x N \9) = Z d N *u e N (p>, s), (27) 

where the nucleon spinor projection is 

2^ u n(p, s)u N {p, s) = — . (28) 

8 Ztj N° 



and the normalization is 



J N 



The nucleon spinor satisfies the Dirac equation with a phase factor associated with the mass 
term due to the CP- violation in the 9 vacuum [21 [3] 

(i 7 • p + m e N e- taNW ' t5 )u e N (p, s) = u e N (p, s)(z 7 • p + m^e"* 01 *^) = 0. (30) 

Now the nucleon matrix element can be written as 

(Ne(p', s')\J* M \N e (p, s)) = u%(p>, s')Vy N (p, s), (31) 

where the vertex can be separated in terms of CP-even and CP-odd form factors 

r! = rr n M) + rfW), (32) 
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Since 9 in QCD is small, one can consider small 9 expansion 
r e ; en (q, 9) = *7^i(<7 2 ) - ^»q»F 2 (q 2 ) + jT e ; en "(q, 0) + ... 

Tf{q, 9) = -<V T5 g / ^ > + 9F' A (q\ 0)(z lfl q 2 - 2m N q,) l5 + -T^'" (q, 0) + ..(33) 
We have used the property that, for small 9 

m% = m N + ^m" N (0)9 2 + ... 

\z e N \ 2 = \z N \ 2 + l -\z N f^)9 2 + ... 

F h2 (q 2 ,9) = F ^{q 2 )+ l -Fl 2 {q 2 )9 2 + ... 
FsAq 2 ,0) = Fi A (q 2 iO)9 + ±F;; A (q 2 ,0)9 3 + ... 



(34) 



The three-point function in Eq. ff26l) is then 



[r7(^) + C(^)l 9F r — • (35) 

N 

Expanding in 9, it gives 

c ™l»-'"M J ), r ,^/e-..''H. e -« { ^i^r r( ,,o)- i " + '"» 



at 2E N 



Ql -Vy-P +m N , odd -1^-p + TJlN 

+ 9[ 2E> T F " M) 2E» 



OCT 1 M W J U J r,7T. OTTl' i M V« > U J T7. J 

+ y G" NJN {q, t f - t > 1, * > 1) + - G^ JJV (g, t/ - < > 1, i > 1) + ...} (36) 

From the current experimental bound on the neutron electric dipole moment, 9 < 10~ 9 
in QCD. Thus, it is sufficient to consider the CP-odd quantities linear in 9. Equating the 
linear 9 term in Eq. (1361) with that in the 9 expansion of the path-integral in Eq. (1251) 



G» JN (q,t,t f ) = (XN{p\t } )J* M {q,t)x N {pMe « (Xn(p, t/)(l + i9Q)jj? M (q } t)x N (p, Q))e=o, 

(37) 

we obtain 

g.mo = o,f,M/-* » i) = |z„| V'.^ g -^{ "' 7 -; F + m " r- fa o) ~ n 9 P F + }, 

(38) 
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and 

G NJN {9 - U,g,t,t/ - t > lj - \Z N \ e Nyt >e { —, 1 (g,U) — 

ia' N (0)m N j 5 , a -z 7 • p + 

-»7-P , + m JV 2 ^(0W 75 

2^ V" ' / 2Ejv / v ) 

a'(0) can be extracted from the nucleon correlator in the 9 vacuum. Consider 

G NN (9,p,t) = (xjv(p,te(p,0)e w «) fl = G NN {0,ftt) + 6G' NN (0,p,t) + 0(9 2 ) + ... (40) 
Taking the appropriate trace and asymptotic limit of the correlator with t 3> 1, we obtain 

Tr(T 4X N(p,t)x N (p,0)) e=0 ee Tr[r 4 G^(0,p,t)^|Z| 2 ^-±^e-^ t 

2777- 

Tr( 75 QxJv(p,t)Xiv(P,0)) e= o ee — iTr [75G jviv (0, p, t) — ^ | Z | 2 — — a' (0) e~ ENt (41) 



A" 

where T 4 = is the projection operator for the time-forward nucleon. If the time t is not 
large enough to filter out the negative parity excited nucleon Su, one will need to use the 
projector 1 + ^p 7 4 where m_/E_ is the mass/energy of Su. 

From the ratio of the two-point functions, one can obtain a'(0), i.e. 

E N + m Tr(~f 5 Qx N (p, t)x N (P, °)) ^qn (A2) 

2m Tr(T 4X N(p,t)xN(P^)) ^ 

In the case of fixed topology that we consider in this manuscript, one can consider the 
two-point functions G V NN in a specific topological sector v 

(43) 

and 



rr[ 75 gW(^ *)]^ \Z\ 2 e- ^ L ~ ^ Xt> ^(0) + ^Tr[ 75 G^(0,p,t » 1)] 



, 1/(1-127/^)^ , -w 

(44) 

With several volumes, one can fit Eqs. (j4"51) and fT44|) and obtain a'(0), \Z\ 2 ,m^, and E^. 

As for the form factors, one can consider the following three-point functions at a fixed 
topology with the sink nucleon momentum p = 0: 

Tr[V 4 G^ j4N (p,tM t>tf ^ 



1 _ u 2/y Xt _ 12 7 /^ X 2 

wxt ±>[[ |lrvMV 



Tr[r 4 G" (0,p,t /; t »!)] + ... (45) 



where Gg(g 2 ) = Fi(q 2 ) + ^2--F 2 (g 2 ) is the electric form factor with q 2 = (En — raw) 2 — <f 
where q = —p. 

x{ie ijk G M (q 2 ) + ~ iU{1 -} 2l,V ^ ) iq % [F*tf, 0) - iAm N F' A {q\ 0) + 2m N F 2 (q 2 )a'(0)}} 
>■ Vxt Am N E N > 

1 — v 2 — Yli IVy 2 

+ 1 ' Xt Tr[iT A >M j G N j ilf ''{0,p,t f ,tZ> 1)] + ... (46) 

*v Xt 

where Gm(? 2 ) = Fi(q 2 ) + F 2 (q 2 ) is the magnetic form factor. 



and 



+ ^ Tr[75r 4 Gjvj4JV (0,p,t/,t » 1)J + ... (47) 

* v Xt 



rr[ ! T l757 ,G^ j4A ,(p,«,i / )] ((/ - ;>>i |Z| 2 e-"»<''- , >e-^' 

x -Ml-127/^?) ^ { j^+^ F , tf;0) + [Fi(?2) + ^m^^,} 
1 _ ^/y^ _ i2 7 /y x g 



+ -^777, J1 ^Tr\ira5jiG' N j 4N (0,p,t f ,t »!)] + ... (48) 



Combining the two-point functions in Eqs. (143]) and (1441) and three-point functions in 
Eqs. (1451) . (j4"6"j) . (|4T|) . and (1451 for several volumes, one can extract the neutron electric 
dipole form factor 

d N (q 2 )=9F^0,q 2 ), (49) 

the anapole form factor 

a N (q 2 ) = 6F^0,q 2 ), (50) 

in addition to the electric and magnetic form factors Ge(g 2 ) and Gm(q 2 ) (and/or Fi(q 2 ) 
and F 2 (q 2 )). When and if the experimental results on the dipole and anapole moments are 
known to be due to the QCD 9 term, one can then determine 9. 

4 Zero Modes, Quenched Approximation, and Anoma- 
lous Magnetic Form Factor 

The role of zero modes in the quark propagators has been discussed extensively in the 
literature. For example, in the quark condensate, its contribution is -^k for a configuration 
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with topological charge v in the 9 vacuum [IT]. Similarly, the pion correlator has a leading 
m 2 V 2 contribution from the zero modes [T8l [T9] . Since one should take the infinite volume 
limit before the chiral limit to have chiral symmetry breaking manifested, the zero mode 
contribution goes away in these limits or their contributions become negligible when the 
volume is sufficient large. However, for CP-odd quantities such as the neutron electric dipole 
moment, the situation is different. We see from Eqs. (1191) . ( H6T) . (H7jL and (HSj) that, for 
CP-odd Green's functions, the leading contribution is proportional to y. Therefore, the zero 
modes will contribute at the finite volume. Specifically, for the case that only one of the 
quark propagators involves zero modes, it would yield a y factor and will, thus, contribute 
to the CP-odd Green's functions to leading order in y. 

For the quenched approximation, the zero mode contribution to the neutron electric 
dipole moment will lead to power divergence in m q at the chiral limit in a finite volume. 
This divergence is not protected by the fermion determinant as is in the dynamical fermion 
case where the determinant measure is proportional to rrv^ Nw and will cancel out the mass 
singularity in the quark propagators due to zero modes. In this case, the divergence is ^3 

for v = 1, and for v > 2 which are inherent in the quenched approximation. Therefore, 
it does not make sense to consider NEDM in the quenched approximation, except perhaps 
to test the algorithm. 



It has been shown [20J recently that, in the light-cone formalism, there is a universal 

relation between the electric dipole form factor and the anomalous magnetic form factor for 
each Fock state, i.e. 



where a denotes the Fock state and (3 a is the P±— and Ti-violating phase in the Fock state a. 
Since /3 a depends on a, F 3 (q 2 ) and F 2 (q 2 ), being the sum of their Fock state components, do 
not have the same q 2 dependence in general, unless the Fock state sum is saturated by a single 
Fock state or [3 a is a pure constant independent of a [21] . We do not see a proportionality 
between F 3 (q 2 ,{)) and F 2 (q 2 ) in the present formulation. Even though the CP-odd three- 
point functions in all v 7^ topological sectors must have the same q 2 dependence in order 
that Eq. ([20} is satisfied, F^(q 2 ) and F 2 (q 2 ) can have different q 2 behaviors in the various 
projected three-point functions considered in Eqs. (H5]) . (pR)j) . (H71) . and (H51) . With precise 
enough lattice simulations, one should be able to check if F 3 (q 2 ) is proportional to F 2 (q 2 ). 
One can also check the conjecture [20J that the neutron and proton electric dipole moments 
repeat the isospin structure of the anomalous magnetic moments, i.e. \d n + d p \ -C \d n — d p \. 

5 Conclusion 

We have studied the Green's function at fixed topological sector for the CP-odd case and 
clarified a question regarding the large V behavior. As an application, we have formulated 
the three-point functions for the electromagnetic current in the nucleon at fixed topology 
which are needed to extract the neutron electric dipole and anapole form factors in this 
practical calculations. It is shown that, in the quenched approximation, the zero modes lead 
to power divergence in the quark mass at the chiral limit for the CP-odd Green functions. 





(51) 
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